Coulomb blockade in a Si channel gated by an Al single-electron transistor 
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We incorporate an Al-AlOj^-Al single-electron transistor as the gate of a narrow (~100 nm) metal-oxide- 
semiconductor field-effect transistor (MOSFET). Near the MOSFET channel conductance threshold, we observe 
oscillations in the conductance associated with Coulomb blockade in the channel, revealing the formation of 
a Si single-electron transistor. Abrupt steps present in sweeps of the Al transistor conductance versus gate 
voltage are correlated with single-electron charging events in the Si transistor, and vice versa. Analysis of these 
correlations using a simple electrostatic model demonstrates that the two single-electron transistor islands are 
closely aligned, with an inter-island capacitance approximately equal to 1/3 of the total capacitance of the Si 
transistor island, indicating that the Si transistor is strongly coupled to the Al transistor. 



Single-electron transistors (SETs) are sensitive electrom- 
eters, able to detect a small fraction of an electron charge, 
with a sensitivity around 10~^ e/VUzi^ They have been ex- 
tensively utilized for metrology^ and may have applications 
in quantum information processing. While metal SETs are 
more common, Si SETs are desirable because of their bet- 
ter stability^ and their ease of incorporation into Si quantum 
computation architectures.^ 

In this letter, we present results from a Si SET self-aligned 
and vertically coupled to a metal SET. An Al-AlO^-Al SET 
is incorporated as the gate of a narrow (~100 nm) ^-channel 
metal-oxide- semiconductor field-effect transistor (MOSFET). 
Near the channel conductance threshold, we observe Coulomb 
blockade oscillations in the conductance, revealing the unin- 
tentional formation of tunnel barriers in the channel and the 
creation of a Si SET. We simultaneously monitor the conduc- 
tance of each SET and observe a clear correlation between 
single-electron charging events in the two SETs, indicating 
that the SET islands are closely situated. The large charge 
(-0.33-0.35 e) induced onto the Al SET island by the addi- 
tion of another electron onto the Si SET island indicates that 
the Si SET is strongly coupled to the Al SET, in the sense 
that the inter-island capacitance is comparable to the total Si 
SET capacitance. Consequently, this device structure could be 
used to study individual tunneling events in the Si SET,^ and 
it represents a single-electron current switch.^ It could also be 
used to probe defect charge motion in MOS structures under 
large electric field at low temperature,^ a potential source of 
decoherence for Si quantum computation.— 

Figure [ita) shows a scanning electron micrograph of a typ- 
ical sample. Fabrication starts with the oxidation of a Si(lOO) 
wafer (resistivity p > 8,000 Q. cm) at 1000 °C, yielding a 
Si02 thickness of about 20 nm. The wafer is selectively ion 
implanted with P at an energy of 50 keV and an areal density 
of 5 X 10 ^^/cm^ to create n-\- contacts. To limit the extent of 
the channel, p-\- regions outside of the n-\- contacts are created 
by another ion implantation of B at 1 8 keV with an areal den- 
sity of 5 X lO^Vcm^ [Fig. [Hb)]. The peak densities of both 
dopants are high enough to conduct at 20 mK and reside close 
to the Si/Si02 interface. If the two implantations are inter- 
changed, a /^-channel device can be made instead, so that both 
polarities can be fabricated on a single chip. After both im- 
plantations, the wafer is annealed at 950 °C for 60 s to acti- 
vate the dopants and to repair implantation damage. Electron- 
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FIG. 1: (a) SEM image of a typical device. The Al SET island 
forms during the second of two evaporations at different angles, (b) 
Schematic of the measurement circuit. The conductance of each SET 
is measured using independent circuits. The red region represents the 
MOSFET conducting channel confined between the two p+ regions. 
The circles containing the letter 'A" represent current- sensitive am- 
plifiers, (c) Coulomb blockade oscillations of the Si SET differential 
conductance as a function of the relative bias y„+ between the Al 



SET and the Si SET at Vo 



-0.604 V. 



beam lithography and self-aligned double angle evaporation 
are used to fabricate the leads and island of the Al SET,^^ as 
well as an Al side gate, used to modulate the conductance of 
both the Al SET and the MOSFET channel [Fig. [Ha)]. Fi- 
nally, the sample is annealed at 425 °C in forming gas for 30 
min to passivate dangling bonds at the Si/Si02 interface. 

All of the measurements that we present here were made 
on a single device at a temperature of about 20 mK. AIT 
magnetic field was applied to keep the Al SET in the normal 
state. The device survived multiple thermal cycles to room 
temperature and displayed only small background charge off- 
set variations between cycles. To avoid confusion, we present 
data from a single cooldown. 

Figure [Hb) shows a schematic of the measurement circuit. 
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FIG. 2: Simultaneously measured conductances of both SETs. [(a) 
and (c)] Coulomb blockade oscillations of the Al and Si SET conduc- 
tances, respectively, at = —0.610 V. [(b) and (d)] Conductance of 
the Al and the Si SET, respectively, vs Vg and y„+ . 

Both n-\- contacts are dc biased at V^+, while the Al SET is 
grounded except for a small dc bias V^s ~ 100 /iV. An ac ex- 
citation Vac = 10 ™s at 46 Hz is applied between the two 
n-\- contacts to measure the channel differential conductance. 
The two p-\- regions are dc biased at potential Vp+ = —0.800 V 
to confine the channel to a small region between them. Fig- 
ure [Ttc) presents typical data for the channel differential con- 
ductance (Gch = ^ch/^ac) as a function of 14+ at constant side 
gate voltage Vg = —0.604 V. The appearance of Coulomb 
blockade oscillations is surprising, because no tunnel barriers 
are deliberately engineered in the channel. 

In order to determine the coupling strength between the 
two SETs and to infer the proximity of the islands, systematic 
sweeps of Vg and Vn^ are performed while the conductances 
of both SETs are measured. Figure [2b) shows the Al SET 
conductance (Gai = /ds/^ds). and Fig.Od) shows the channel 
SET differential conductance versus Vg and V„+. If each SET 
were electrostatically isolated from the other, its conductance 
maxima would trace out straight lines in these graphs. Discon- 
tinuities from this straight-line behavior evident in Figs. [2b) 
and [2d) are a signature of discrete charging events close to 
the SET islands, events we would expect to observe if the two 
islands were in close proximity. To confirm this hypoth- 
esis, the maxima in Figs. [2b) and [2d) are fitted with Gaus- 
sians, and the resulting peak centroids are plotted in Fig. [3] 
The two SETs display a clear correlation: whenever an Al 
SET conductance peak trace meets one from the Si SET, it 
makes an abrupt step, and vice versa. Because each peak trace 
corresponds to a unit change in the number of electrons on the 
corresponding SET island, this correlation proves that single- 
electron charging events in one SET are coupled to the other. 

To more quantitatively explain the above results, we have 
modeled the device using the circuit depicted in Fig. [Ua). 



Under the assumption that the system can minimize its elec- 
trostatic energy automatically by independently adjusting the 
number of electrons A^si and A^ai on the two SET islands, the 
total electrostatic energy of this circuit is given in matrix form 
byi4 



E{NAuNsi,Vg,Vn^) = -Q 



2- V 



G. (1) 



Here, Cz_a\ = Ci + C3 + Q + Q_ai and Ci:_si = C2 + C4 + 
Cc + Cg_si are the total capacitance of the Al SET and of 
the Si SET island, respectively, and Q = {—cNai + Q_Ai^g + 
C^Vn^^—eNsi + Cg_siVg + C2Vn^y includes the virtual and ac- 
tual charges on the SET islands. Under energy degenerate 
conditions. Coulomb blockade is lifted, resulting in the max- 
imal SET conductances. There are in total six such degener- 
acy conditions associated with adding or subtracting one elec- 
tron from an SET island, determined by E {Nai , Nsi , Vg , y„+ ) = 
E{NAi^5NAuNsi^5NsuVg,Vn^)^{5Nsi)eVn^. Here, 5A^ai = 
0,±1; 5A^si = 0,±1; and |5A^Ai+5A^si| < 2, and (5A^si)^K+ is 
the extra work done by voltage source y„+ when one electron 
tunnels through junction €2-^^ These equations establish the 
hexagonal phase diagram depicted in Fig. [Jb), and there is 
a correspondence evident between this diagram and the data 
in Fig. [3l However, the capacitances associated with the Si 
SET appear to be bias voltage dependent, resulting in the non- 
identical hexagons in the data. 

There are in total seven capacitance parameters in our cir- 
cuit model. Ci = 282 ± 6 aF is extracted from diamond chart 
measurements (not shown) with the MOSFET in the off state 
(Vn^ = V). Ci is dominated by overlap between the Al SET 
leads and its island and should be insensitive to the pres- 
ence or absence of an underlying MOSFET channel. The 
remaining six parameters can be extracted from the slopes 
5"!, 5*2, S3 and the separations Ai, A2, A3 of each hexagon 
in Fig. [3 as defined in Fig. ^h). We extract these param- 
eters from the data as follows. First, the boundaries given 
by SNai-^SNsi = ±1 (the nominally straight lines traced out 
by the data in Fig. [3]) are each fitted to a line. Boundaries 
corresponding to 5A^Ai+5A^si = (an effective transfer of an 
electron from one island to the other) are not clearly visible, so 
they are determined by neighboring intersections of the visible 
boundaries. To compensate for gradual changes in the capaci- 
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FIG. 3: Conductance maxima of both SETs vs y„+ and Vg. Red 
dots and blue dots are Gaussian fits to the data in Figs. [2b) and 
[2d), respectively. Black lines are a linear fit to the points on each 
edge. The regions labeled a, b, c, and d are the four hexagons whose 
parameters are presented in Table [H 
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FIG. 4: (a) Circuit model for the coupled SET system. A^ai ^nd A^si 
are the number of electrons on the Al and the Si SET island, respec- 
tively. Due to the very small drain- source bias of each SET, we can 
simplify the two tunnel barrier capacitances for each SET to a sin- 
gle capacitance (Ci and C2) as shown, (b) Hexagonal phase diagram 
based on the model in (a). Each hexagon represents a configuration 
with a different number of charges on the SET islands. Si, S2, and 
^3 are the slopes of the hexagon edges. Ai, A2, and A3 are the sepa- 
rations between opposite parallel edges of the hexagon. 



TABLE I: Capacitances of the four hexagons labeled in Fig.[3]for the 
circuit model in Fig.lH 



hexagon 


C2(aF) 


Q(aF) 


C3(aF) 


C4(aF) 


C,_Ai(aF) 




a 


20±6 


16±1 


I4±l 


7±1 


4.1±0.2 


6.0±0.2 


b 


32±6 


21±1 


7±1 


1±1 


4.6±0.3 


6.0±0.3 


c 


32±4 


22±1 


7±1 


2±1 


4.1±0.3 


6.7±0.2 


d 


31±5 


21±1 


10±1 


2±1 


3.9±0.2 


5.9±0.2 



tances with bias voltage, averages are made for the slopes and 
separations from opposite boundaries within each hexagon. 
Then, the six unknown capacitances can be extracted by solv- 
ing six analytical equations relating 5"!, 5*2, S3, Ai, A2, and A3. 

Discrete background charge motion near the SET islands, 
which changes the electrostatics of the system, makes sys- 
tematic study of all the hexagons in Fig. [3] difficult. The 
capacitances for the four typical hexagons labeled in Fig. [3] 
are presented in Table HI For hexagon a, C^.si = 49 ± 6 aF. 
A diamond chart measurement of the Si SET near this bias 
point (not shown) gives C^.si = 49 ± 3 aF, in good agree- 



ment and confirming the validity of the circuit model. Signif- 
icantly, for all the hexagons in TableU Q/Ci:_si ^ 33%-35% 
indicates that the Si SET is strongly coupled to the Al SET, 
while Cc/Ci^_A\ ~ 5%-7% explains why the discontinuities in 
Fig. [2d) are less obvious than those in Fig. [2b). 

If the overlap between the two SET islands were perfect, 
the value of Q as calculated from the Al SET island dimen- 
sions and the Si02 thickness would be about 30 aF. This is 
close to the values in Table [J for hexagons b-d. The small 
values for C4 in these hexagons mean that there is almost no 
overlap between the Si SET island and the Al SET leads. This 
strongly suggests that the induced Si SET island is located 
directly beneath the Al SET island. 

We hypothesize the following reasons for the formation of 
an aligned SET in the channel. Although the width of the SET 
island and that of the leads in the evaporation shadow mask are 
the same, the SET island is formed during the second evapora- 
tion. A slow pinch-off of features in the mask during the first 
evaporation therefore makes the island slightly narrower than 
the leads. If the angle between the evaporations is incorrect, 
there may also be a lateral offset between the island and leads. 
This island/leads width asymmetry and lateral offset may lead 
to lateral constrictions in the MOSFET channel below, creat- 
ing tunnel barriers and therefore an SET in the channel aligned 
with the Al SET above. 

This SET sandwich architecture could be used to char- 
acterize the MOS structure at low temperature via a cross- 
correlation measurement between the two SETs.^ For exam- 
ple, measurements of the Al SET show that about 10 electrons 
have already accumulated on the Si SET island when the first 
measurable Si SET conductance peak appears. This architec- 
ture could help identify sources of unwanted charge motion 
that may also be sources of decoherence for Si quantum com- 
putation. Because the Si02 layer could be made much thinner, 
future experiments could more fully explore the strongly cou- 
pled two- SET regime. 
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